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Abstract 

Thermal instability is examined for advection-dominated one-temperature 
accretion disks. We consider axisymmetric perturbations with short wavelength 
in the radial direction. The viscosity is assumed to be sufficiently small for the 
vertical hydrostatic balance to hold in perturbed states. The type of viscosity 
is given either by the a-viscosity or by a diffusion-type stress tensor. Optically 
thick disks are found to be in general more unstable than optically thin ones. 
^ When the thermal diffusion is present, the optically thin disks become stable, 

00 but the optically thick disks are still unstable. The instability of the advection- 

^ dominated disks is different from that of the geometrically thin disks without 

^ advection. In the case of no advection, the thermal mode behaves under no 

0\ appreciable surface density change. In the case of advection-dominated disks, 

however, the thermal mode occurs with no appreciable pressure change (compared 
I with the density change), when local perturbations are considered. The variations 

of angular momentum and of surface density associated with the perturbations 
^ lead to a thermal instability. The astrophysical implications of this instability 

are briefly discussed. 
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1. Introduction 

Thermal instability of geometrically thin accretion disks has been extensively studied 
with fruitful applications to many time varying phenomena in dwarf novae. X-ray stars, etc. 
The mechanism of this thermal instability is well understood, and its essential part can be 
explained by considering perturbations which occur without changing the surface density 
of the disks (e.g., Pringle 1976). The fact that this is a good approximation is closely 
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related to the fact that the disk is geometricaUy thin. In geometricaUy thin disks, the 
radial force balance is realized by the balance between the centrifugal and gravitational 
forces. Because both forces are much stronger than the pressure force, a temperature 
change associated with the thermal mode occurs with no motion in the horizontal plane. 
The gas expands or shrinks only in the vertical direction so that a hydrostatic balance in 
the vertical direction is realized without changing the surface density. This approximation 
is valid even when perturbations are local in the radial direction. 

In advection-dominated disks, however, the situation is different. Optically thick disks 
of this type have been introduced by Abramowicz, Lasota, & Xu (1986) and studied in 
details by Abramowicz et al. (1988), Kato, Honma, & Matsumoto (1988), and Chen & 
Taam (1993). Optically thin advection-dominated disks have been most recently inves- 
tigated independently by Abramowicz et al. (1995) and Narayan & Yi (1994, 1995a, b) 
and by Chen et al. (1995). [It is noted, however, that optically thin advection-dominated 
transonic disks are already introduced by Matsumoto, Kato, & Fukue (1985), although 
full details are not published.] Advection-dominated disks are not vertically thin since 
the local disk half-thickness H(r) is close to the radial coordinate r itself (Abramowicz et 
al. 1988, 1995; Narayan & Yi 1995a). This means that the pressure force in the radial 
direction is non-negligible compared with the centrifugal and gravitational forces. There- 
fore a pressure change associated with the thermal perturbation brings about a motion 
in the radial direction. This results a variation of surface density which is larger than 
the pressure variation when a local perturbation is considered, as will be discussed later. 
Thus, the approximation of no change of surface density is no longer valid and the criterion 
of thermal instability becomes quite different from that in the case of geometrically thin 
disks, especially when the local perturbation is considered. This difference has not been 
fully recognized in the previous studies of thermal instability of advection-dominated disks 
(Wallinder 1990; Honma et al. 1991; Chen & Taam 1993). A purpose of this paper is 
to clear this difference and to derive the instability criterion against local perturbations. 
Local perturbations mean here that the radial wavelength of perturbations is shorter than 
r, where r > H . 

In addition to the above formal interest in thermal instability of advection-dominated 
disks, this kind of instability may have important astrophysical applications. It was 
recently recognized that advection-dominated disks are possible models of some AGN 
(Abramowicz et al 1995), of X-ray transients (Narayan & Yi 1995), and of the Sgr A* 
(Narayan, Yi, & Mahadavan 1995). Some of the time variations observed in these objects 
will be related to the thermal instability of advection-dominated disks. 

2. Vertically Integrated Equations 
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We consider axisymmetric perturbations on axisymmetric disks. To this examination 
we introduce a cylindrical system of coordinates (r, (^,2;) which is centered on the central 
object and its 2;-axis is chosen as the rotation axis of the disk. 

2.1. Validity of Vertically Integrated Equations 

We focus here the advection-dominated accretion disks. Considering that H/r is of 
order of unity and also the presence of the fast radial flow motion, we must be careful 
about the validity of using the vertically integrated equations which are widely used in the 
analyses of thermal instability of geometrically thin disks. 

To examine the validity condition we compare here the time scale associated with 
perturbations with the dynamical time scale in which the vertical hydrostatic balance is 
realized. The growth rate of the thermal instability in geometrically thin a-disks is of the 
order of aO, where O is the angular velocity of disk rotation. However, as we show later, in 
advection-dominated disks, for perturbations with short radial wavelengthes, the growth 
rate becomes aQA. The factor A depends on the model of the dissipation processes. In 
the case which will be discussed in section 4, we have A = (kry/'^(H/r)^ where A = 1 and 
A = {kHY in the cases discussed in sections 5.1 and 5.2, respectively. Here k is the radial 
wavelength of perturbations. In addition to this, we must consider the time variation due 
to the propagation of the perturbation. The frequency associated with this is kU ^ where 
U(r) is the radial velocity of advection flow, 

V~.ar(^^)\ (1) 

and Cg is the local sound speed. These frequencies must be lower than the dynamical 
frequency, Cs/H(= O), in order that one may use the vertically integrated equations since 
they are derived by assuming hydrostatic balance. The condition is 

aOA < O and kU < O, (2) 

and can be rewritten as 

a < A~^ and a<(kH)~^^. (3) 

H 

Inequality (3) shows that local perturbations, A;r ^ 1, can be studied using the verti- 
cally integrated equations as long as a is sufficiently small. In this paper we restrict our 
discussion to this case. 

2.2. Basic Equations 
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We assume the disk is axisymmetric and non-self-gravitating. The basic equations 
which we adopt here are the time-dependent version of those given by Matsumoto et ah 
(1984). That is, the basic hydrodynamic equations are integrated in the vertical direction 
under the assumption that a hydrostatic balance always holds in the vertical direction. To 
avoid ambiguities regarding the vertical integration, the integration is carried out under 
the assumption that the pressure, p, is related to the density, /?, through the polytropic 
relation, i.e., p oc p^~^^/^ ^ with iV = 3 for simplicity. The effects of general relativity are 
introduced by the pseudo-Newtonian potential, tp = —GM/ (R — r^), of Paczyhski & Wiita 
(1980). Here R is the distance from the central object of mass M and Vg = 2GM I is its 
Schwarzschild radius. Under these assumptions, the equation of continuity is 

f + i|:(^«0 = o. ,4, 

where E is the surface density obtained by the vertical integration of /?, and U{< 0) is the 
radial component of the velocity. Hereafter, U is taken to be negative, although U used 
before was the absolute value of the radial flow. The r- and components of equation of 
motions are, respectively, 

DU ^2 I dW W c/IuOk 



DV UV Id 



Dt r ^ dr ^ dr 

and 

; ) \/ I I \/ 1 'ri ^ 



_ + + -^^{r'Wr^) = 0, (6) 

Dt r r^T, or 

where W is the vertically integrated p, V is the (^-component of the velocity, and OK(r) is 
the Keplerian angular velocity on the equatorial plane, i.e.. 



K 



1 dt/j 
r dr 



1/2 



(7) 

z=0 



The operator D / Dt is deflned as 



In the case of iV = 3, the vertical integration gives (Hoshi 1977) 

E = 2hp,H, W = 2hp,H, (9) 

where H is the half disk-thickness, pc and pc are the density and the pressure on the 
equatorial plane respectively, and Is = 16/35 and I4 = 128/315 are numerical constants. 
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The quantity W^p in equation (6) is the vertical integral of the r(^-component of the viscous 
stress tensor. We adopt here the standard a-viscosity (Shakura & Sunyaev 1973) 

Wr^ = aW. (10) 



The hydrostatic balance in the 2;-direction gives 

{QkHY =2{N + 1)^. (11) 

pc 

In the case of iV = 3, the vertical integration of energy equation can be written in a simple 
form (Matsumoto et al. 1984): 

|eT.^ = Q+-Q-. (12) 

Here, Tc and Sc are the temperature and the specific entropy of the gas on the equatorial 
plane respectively, is the heating rate per unit surface given by 

Q+ = -rlF.,|:(^), (13) 

and Q~ is the radiative flux escaping from the disk surface. The detailed expression for 
Q~ will not effect the advection-dominated disks and thus is omitted here. Finally, the 
equation of state for optically thin disks is 

Pc = —^PcTc. (14) 

For optically thick disks the radiation pressure, (l/3)aT;,'^, is added on the right-hand side 
of equation (14). 

In section 5, some modiflcations of these equations are discussed: Instead of the 
standard a-viscosity, a diffusion-type stress tensor is adopted. Furthermore, i) a viscous 
force is added to the right-hand side of equation (5) and ii) a thermal diffusion term is 
taken into account on the right-hand side of equation (12). 

3. Equations of Perturbations under Geostrophic Approximation 

A small-amplitude, short-wavelength perturbation is superimposed onto an unper- 
turbed disk. At the radius where local perturbations are imposed, the unperturbed disk 
is assumed to be rotating cylindrically, O = 0(r), and the radial velocity of the flow Uq 
depends only on the radius, Uq = Uo(r). Eulerian perturbations superimposed over the 
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unperturbed quantities Eq, Uq^ Vo(= fir), and Wq are denoted by Ei, ?7i, Vi and VFi, 
respectively. Hereafter we adopt dimensionless variables: 



Si W^i ^1 , ^1 ...^ 

(7 = —, ti^ = T7r' ^ = TT "=rr- 
Zjq Wo Uo ilr 

They are taken to be proportional to exp(nOt — ikr)^ where n is the dimensionless growth 
rate of perturbations normalized by O. Hereafter, the subscripts to E, VF, and U are 
omitted for simplicity. 

From equation (4) we have 

n*cr — ikr — u = 0, (16) 
ilr 

where ^ 

= n — ikr——. (17) 
ilr 

When considering the r-component of equation of motion, we introduce an approximation 
to filter out the acoustic waves. For low frequency perturbations in a rapidly rotating 
system, the pressure force is roughly balanced by the Coriolis force. That is, the term 
of DU I Dt in equation (5) can be neglected. This approximate procedure is well known 
as the geostrophic approximation in meteorology and oceanography (e.g., Pedlosky 1979). 
(We can show later, using the final results, that the terms neglected are higher order small 
quantities.) Under this approximation, we have from equation (5) 

( c^V ^ln(iyOK) /cs V 
-^v-^kr[^-^ w - ^^^^ J a = 0, (18) 

where Cg is the speed of sound defined by = VF/E. The (^-component of equation of 
motion (6) leads to 



-*^ + ^o;^" + "^-^"^'(vo;J " = °' ^^^^ 

where k is the epicyclic frequency defined by = 20(20 + rdQ/dr). Finally, energy 
equation (12) yields, after lengthly calculations, 

n^,Q(Ciw — €2(7) = GvV + GwW + G„a + (20) 

where 

Gv = iaQkr (21a) 
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Gu = {21d) 

Here, and Q~ are the viscous heating and radiative coohng rates in the unperturbed 
state, and is the advection coohng rate in the unperturbed state, i.e., 

Q:, = |ETe^ = Q+-Q-. (22) 

The first term on the right-hand side of equation (20) comes from variation of heating due 
to a change of rotational velocity V [see equation(13)]. In the above equations C's are 
numerical constants defined by 

2(r3-i)' 2(r3-i)' ^ ^ 

and 

- 2(4-3/3)' - 2(4-3/3)' ^^^^^ 

where Fi and Ts are adiabatic indices which include the effects of radiation pressure, and 
/3 is the ratio of the gas pressure to the total pressure. 

In the case of optically thick, advection-dominated disks, we have 7 = 4/3 and /3 = 0, 
since pr ^ Pg- Then the constants defined above become 

Ci = l, C^ = \ C^ = \, C, = \. (24a) 

On the other hand, in optically thin, advection-dominated disks, we have 7 = 5/3 and 
/3 = 1, since the disks are supported by ion pressure. Hence we adopt 

Ci = 2, C2 = 3, C3 = 1, C4 = 0. (246) 

4. Thermal Instability agaist Local Perturbations: a viscosity 

Before examining in detail the thermal instability in advection-dominated disks, we 
shall discuss briefiy the reason for thermal instability in geometrically thin disks. Compari- 
son between these two is helpful to understand the essential difference of thermal instability 
in geometrically thick disks and geometrically thin disks. 
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4.1. Local Instability of Geometrically Thin Disks 



It is well known that the thermal instability in geometrically thin disks can be exam- 
ined under the approximation that there is no surface density change during the growth of 
the instability. Under this approximation, the v and a terms in the energy equation (20) 
can be neglected. Furthermore, the GuU term in that equation can also be neglected, since 
Gu is negligibly small in disks with no advective energy transport. Thus, energy equation 
(20) can be approximated as 

nQCiw = GwW, (25) 

where 

Here, term has been neglected in G^, since the advective cooling is negligible in 

geometrically thin disks. The condition of the thermal instability is then 

G^ > 0. (27) 

This is the well-known criterion for thermal instability (e.g., Pringle 1976). Equation (25) 
shows that the order of n is a. 

If the a terms are neglected in the r- and (^-components of momentum equation, we 
have from equations (18) and (19) 

-2v - ikr(^^ u; = 0, (28) 
\\lr J 

and 



+ TTTTT^^^ ~ iakri -— ] w = 0. (29) 
202 Or J V ; 

These equations show that 0(v) = (kH)(H /r)O(w) and 0(u) = (kr)O(w)^ where O de- 
notes the order of the magnitude of the subsequent argument. Furthermore, the continuity 
equation (16), 

ncr — ikr——u = 0, (30) 
gives 0(a) = {kH){H j r)0{u). Thus, with kH and Hjr < 1, we have 

0{v) < 0{a) < 0{w) < 0{u). (31) 

In this way, the neglection of the cr-terms in the momentum equation is found to be consis- 
tent with the final results. The above results also confirm the validity of approximations 
used to derive equation (25). 
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After obtaining n from equation (25), we can derive u, and cr, successively, from 
equations (28) - (30), as functions of n and w. We shall discuss now how the above 
situations are changed when local perturbations (kr ^ 1) are considered in advection- 
dominated disks (H < r). 

4-2. Local Perturbations in Advection- Dominated Disks 

Advection-dominated disks are not geometrically thin, H/r < 1. This means that 
the pressure force is non-negligible compared with the gravitational and centrifugal forces. 
Hence, if a temperature increment occurs at a certain localized region, the disk at the 
region expands in the radial direction. For a fixed magnitude of radial expansion, the 
surface-density decrease due to the expansion is larger in the case when radial wavelength 
is shorter. On the other hand, the decrease of vertically integrated pressure in that region 
due to the radial expansion is smaller when the radial wavelength is shorter. This follows 
from the fact that the pressure at a particular radius is determined by the weight of the mass 
supported by the effective gravitational force (the difference between the centrifugal force 
and the gravitational one) outside this radius. (This can be understood by integrating the 
r-component of equation of motion from r = oo to the radius in question.) In other words, 
the pressure can not change much (compared with the density change) by local perturbation 
since it does not change the global structure of the force balance. In summary, we can 
expect a large density change, but a small pressure change, when local perturbations are 
considered. That is, 0(w) < 0(a). This is a difference from the case of geometrically thin 
disks. 

The fact that in local perturbations a large density change occurs with a small pres- 
sure change is well known in stellar hydrodynamics in relation to flickers. A flicker is a 
phenomenon which occurs by thermal instability of a localized perturbation in shell burn- 
ing region of stars, found flrst by Schwarzschild & Harm (1965) (e.g., Fujimoto, Sugimoto 
1979). The thermal instability of geometrically thick disks agaist local perturbations is 
similar to the flickers in the sense that 0(w) < 0(cr), although the heating sources leading 
to instabilities are different in each case. 

Let us now consider the orders of variables, u, and a. In the present case the 

order of is a(kry/'^(H/r)^ as will be conflrmed later. If this is adopted, the equation 
of continuity, equation (16), shows that the amplitude of u is smaller than that of a by 
a factor of (kr)~^/'^(H/r)~^ ^ since the order of U is given by equation (1). That is, we 
have 0{u) = {kr)-^/^{H/r)-^0{a). This means that the u term in equation (19) can 
be neglected, compared with the term of a. Comparison of equation (18) and equation 
(19) shows further that the flrst term on the left-hand side of equation (18), — 2u, can be 
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neglected. [The w and a terms in equation (18) have the orders of the corresponding terms 
in equation (19) times . On the other hand, the v term in equation (18) has the order 
of the corresponding term in equation (19) times a~^{kr)~^/'^{H /r)~^ ^ which is smaUer 
than Thus, equations (18) and (19) can be approximated, respectively, as 



-'^'[sTr) " - dlnr [sTr ) " = °' ^''^ 
and 

-*^ + ^o;"-^"^'(,o;J " = °- ^^^^ 

These are the r- and (p- components of equation of motion. 

We have mentioned that in local perturbations in advection-dominated disks the order 
of w is smaller than that of a. Equation (32) really shows that the ratio is about (kr)~^ , 
i.e., 0(w) = (kr)-^0((7). Equation (33) shows then that 0(v) = (kr)-^/^(H/r)0((7). In 
summary, we have 

0{w) < 0{v) ~ 0{u) < 0{a). (34) 

The use of the above relations among u, and a leads us to an approximate form of 

equation (20): 

n*0( — (72cr) = GvV. (35) 

This is the energy equation. Combination of equations (32), (33) and (35) gives the final 
dispersion relation. 



2 a Gi 



fi 

* C2 o \ nr 



dln{WQK) ^ «^ U 



dlnr 202 a0r(cs/0r)2 

dln{WQK) U 



+ 



(36) 



dlnr 202 a0r(cs/0r)2 

This shows the presence of an unstable mode whose growth rate is of the order of aO(A;r)"^/2 
(H/r). This growth rate is consistent with the initial assumption concerning the order of 



The right-hand side of equation (35) shows that the direct source of the thermal insta- 
bility is angular momentum variation. The angular momentum variation occurs as a nat- 
ural consequence of density variation. If angular momentum per unit surface is conserved 
(in the Lagrangian sense) during the growth of perturbations, we have [dCEV) / dt]i = 0. 
This equation relates directly v and a when the change of ?7, i.e., u, is neglected. That is, 
we have -|- (k,'^ /2Q'^)(U/Qr)a = 0. This is an approximate version of equation (33). If 
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this equation and equation (35) are combined to derive a dispersion relation, we obtain an 
equation which is the same as equation (36) except that the first term in the brackets of 
equation (36) is omitted. 

^.3. Comments on the Case with no Shear 

As shown above, the main source of the thermal instability in geometrically thick 
a-disks is the angular momentum variation associated with perturbations. However, con- 
sideration of the case when the effect of angular momentum variation is neglected seems 
to be instructive for understanding the relation of the present instability with that of non- 
rotating stars (i.e., flickers). If the effects of angular momentum variation is neglected, the 
main term on the right-hand side of equation (20) is only Ga-c since G^'w and GuU can 
be neglected simply because 0{w),0{u) < 0(a). Then, one can obtain easily the growth 
rate as 

G 

n*0 = -7^, (37) 
and the condition of instability is thus Go- < 0. 

This instability criterion can be understood as following. Let us consider a surface 
density decrease at a localized region. If Go- < (which is true in most of the cases in 
advection-dominated disks, see section 5), this density decrease leads to an entropy increase 
at the place, since the right-hand side of equation (20) represents entropy change. Since 
an entropy change is proportional to Ciw — 02(7 ~ — (72cr, the entropy increase results in a 
decrease of the surface density. Therefore, the surface density decreases further more. This 
instability mechanism is the same, in principle, as that of flickers in stars. (It should also 
be noted that in the present case the term of —2v should be retained on the left-hand side 
of equation (32), but this makes no essential difference in the argument presented above.) 

5. Thermal Instability against Local Perturbations: Diffusion- Type Viscosity 

5.1. The Case with no Thermal Diffusion 

The arguments in the previous section have some drawbacks. First, the results show 
that the spatial variations of angular momentum due to the perturbation are the main 
source for thermal instability. This suggests that we must re-examine the same problem in 
the case where the stress tensor has a diffusion form, because the stress tensor described 
by the a-model does not reflect properly the effects of variations of angular momentum. 
As the r(^-component of the stress tensor we adopt here 

'"- = -'"■*: (7)' 
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where rj is the verticaUy integrated turbulent dynamical viscosity. Its magnitude is about 
a*SCg/0 and a* is the dimensionless viscosity parameter corresponding to the a used 
before. If equation (38) is used instead of equation (10), we must change equation (19) to 



^2 jj 



202 Q 



r 



\lr 



c/lnO / dlnrj 



d\nr V c?lnVF 



w + 



c/lnO / dlnrj \ 



dlnr 



V ainE ) 



a 



tftrv 



w 



(39) 



Second, in the conventional a-model the viscous force in the radial direction is not 
taken into account in the r-component of momentum equation. However, the force might 
have non-negligible effects on behavior of thermal perturbations, since we are considering 
perturbations whose wavelength in the radial direction is short. The main term which 
should be added on the right-hand side of equation (5) is rjd^U / dr^ . Hence, equation (18) 
is changed to 



-2v 



w 



dlnr V Or 



a 



u. 



(40) 



In addition to the above modifications, equation (20) must also be changed to 

n^,Q(Ciw — €2(7) = GvV -|- GwW -|- G„a -|- (41) 



where 



Q+ f d\n'q\ Q-fd\nQ-\ ^ Q'^ 



W [dlnWj^ W [ dlnW )^ if ' ^^^^^ 

/ ain, \ Q- fdlnQ- \ Q^, 

Gu = -^- (42rf) 

The last basic equation, equation (16), remains unchanged. The set of equations to be 
used here is equations (39) - (41), and (16). 

The main terms in equation (39) are those in the brackets, since they are multiplied 
by the factor kr. This means that the angular momentum transport by viscosity is so 
efficient that d{r'^Wircp) / r'^ dr ~ holds. Furthermore, among the terms in the brackets. 
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the last two are the dominant terms. Thus 0{v) = (kr)~^ 0(a). The final results which 
will be given below show that in the present case the order of is a*, i.e., 0(n*) = a*. 
When this is taken into account, we have 0{u) = (kr)~^ (H / r)~'^ 0(a) from equations (16) 
and 0{w) < (kr)~^ 0(a) from equation (40). In summary, we have 

0(v) ~ 0(w) ~ 0(u) < 0(a). (43) 

Hence, the basic equations (16), (40), (39), and (41), describing the local thermal pertur- 
bations, can be approximated, respectively, as 

n*cr — ikr——u = 0, (16) 
Or 

dln(WnK) ,2 

—ikrw a = — a*( fcrj u, (44) 

amr 

c/lnO / dlnr] \ ^ 

and 

n^Q(-C2a) = G^v + G^a. (46) 



The elimination of v from equations (45) and (46) gives 

nM-C2(j) = Ga (47) 

where 



- 1 (/In0/aim7^\ - 
Q+ ( dlnr, \ _Q2fdlnCr\ 



ad (4-9,) 

and the condition of instability is 

G <0. (49) 

The growth rate is given by —G/C2 and which is of the order of in consistent with the 
assumption made before. Equations (16) and (44) are subsidiary relations which express 
u and w in terms of a and n*. 

Among the three terms on the right-hand side of equation (48), those with and 
are the main terms, since ~ ^ Q~ in advection-dominated disks. In 

the conventional model of turbulent viscosity, we take r] oc Ec^/O oc W/Q and thus 
(dlnrj / dln'S)w ~ 0. This means that the term of is small and G ~ G^- Hence, 
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the growth rate is practicaUy equivalent to that given by equation (37). In opticaUy thick, 
advection-dominated disks, is negative since > 0. In the case of opticaUy thin, 
advection-dominated disks, however, the major term of Go- vanishes in the lowest approx- 
imations since C4 ~ 0. This means that in the case of optically thin disks, some of the 
terms neglected in deriving equation (48) should be restored to derive the growth rate. 
However, the main effect of the inclusion of these terms is just adding an imaginary part 
to n*. The real part of is still given by equation (48), unless G is too small. Hence, the 
instability criterion is still given by G < 0(or Go- < 0), and the disks are unstable, since 
(dlnQ~ / dln'E)w > for the free-free cooling. 



5.2. Effects of Thermal Diffusion 

In the arguments made in our discussion so far, an important process was still omitted. 
Since we are interested in perturbations whose radial wavelengths are short, the thermal 
diffusion in the radial direction cannot be neglected. That is, on the right-hand side 
of energy equation (12), the term of vertically integrated div(A'gradT) should be added, 
where K is the thermal conductivity and T the temperature. Since turbulence is considered 
here as the origin of viscosity, K/C^ (Cp is the specific heat of constant pressure) must 
be at least of the order of the eddy viscosity, i.e., K/Cp ~ a^pCsH ~ a^QpH^. Hence, 
the vertical integration of K is of the order of a*0Eif2(VF/ET) = a^QH^W/T. When 
perturbations with a short wavelength in the radial direction are considered, the perturbed 
part of the vertical integration of div(A'gradT) is of the order of — a*0(A;if)^VF(Ti/To)o, 
where Ti is the temperature perturbation over the unperturbed one To , and the subscript 
outside the parentheses means the values at the equator. Hereafter we write the quantity 
as — a*0(A;if)^VF(Ti/To)o, where a* is of the order of a*. If the vertical hydrostatic 
equilibrium is assumed even in the perturbed state, we obtain from the variation of equation 
of state 

(^) = + (Ci - l)a. (50) 

The above consideration shows that in the present case equation (41) can still keep 
its form, if G^ and Go- are modified as 
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The last terms of equations (51a) and (51b) are added to equations (42b) and (42c), 
respectively. Local perturbations in geometrically thick disks imply {kHY ^ 1, and 
therefore the added terms are the main terms of and Go-- 

Again, equations (39) and (40) suggest 0{v) = (kr)~^ 0(a) and 0{w) = (kr)~^ 0(a). 
Also, 0(n*) = d^,(kH)^ implied from the final results and equation (16) leads to 0{u) = 
(kr)O(a). In summary we have 

0{w) ~ 0{v) < 0{a) < 0{u). (52) 

On the right-hand side of energy equation (41), Ga is larger than Gu by a factor (kH)^ . 
This means that in equation (41) GuU can be neglected in comparison with Ga-a^ although 
u is larger than a. A similar argument can be made to neglect the terms of GjjV and G^w. 
These considerations lead us to the following approximate set of equations: 

n^a — ikr — u = 0, (16) 
ilr 

dlniW^K) ,2 

—ikrw a = — a*( fcrj u, (44) 

amr 

U , / Cs \ ^ c/lnO / d\nr]\ , , 

Yn,^.,mn.^^_^u^...Mr[-^) _(^_|j^^„ = 0, (53) 

and 

nM-C2o) = G^a. (54) 



Equation (54) shows that the condition of instability is 

G^ < 0, (55) 

and the growth rate is given by — Gct/C*2, where Go- is given by equation (51b), not equa- 
tion (42c). The remaining equations (16), (44), and (53), are subsidiary relations which 
determine u, and w for given a and n*. Inequality (55) implies that the condition of 
instability is practically C4 > 1, since (kH)^ > 1. This condition is equivalent to 

f^<l- (56) 

Thus, optically thick advection-dominated disks (which have /3 0) are unstable, while 
optically thin ones (which have /3 1) are stable in the limit of short wavelength pertur- 
bations. 
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The cause of this instabihty can be understood as foUows. Let us assume a surface 
density decrease at a certain radius. If 6*4 > 1, that leads to a temperature decrease 
[see equation (50)]. Then, heat flows to the place by thermal diffusion process, leading to 
an increase of entropy there. Since Si oc Ciw — (72cr (where Si is the entropy change), 
an entropy increase brings about a decrease of the surface density. This acts to amplify 
the initial surface density decrease and sets up an instability. The above argument shows 
that the essential point of the instability is a positive correlation between the changes 
of temperature and surface density, i.e., C4 > 1. This positive correlation is realized 
in the radiation pressure dominated disks. In that case, since the hydrostatic balance 
in the vertical direction implies p/p oc VF/E oc if^, we have T oc p^^'^ oc (Eif)"^/"^ oc 
[S(VF/E)"^/^]"^/'^ oc VF-^/^E"^/^. So the changes of temperature and surface density occur 
in the same direction, implying C4 > 1. On the other hand, in the case of no radiation 
pressure we have T oc p/ p oc VF/S, and the variations of T and E occur with the opposite 
sign, i.e., < 1. 

6. Discussion 

We have discussed the thermal instability of advection-dominated disks against local 
perturbations. In the case a-viscosity, the disks are always unstable [see equation (36)], 
while in the case a diffusion-type viscosity they are unstable if a certain condition is 
satisfled [see equation (49) or equation (55)]. The essence of this instability is related to 
the fact that such disks have a non-negligible pressure force compared with the centrifugal 
and gravitational forces. Because of the pressure force, the gas can expand (or contract) 
in the radial direction when a perturbation is imposed. If the radial wavelength of the 
perturbation is short, we can expect a large density change by the radial expansion (or 
contraction). This can occur without much pressure change. That a large density change 
associated with a small pressure change is the essential ingredient of the instability. The 
energy sources giving rise to the instability are the unbalance of the heating and the 
cooling rates caused by a density change and a spatial shear variation associated with the 
perturbation. 

We should emphasize here that we have examined the thermal instability of disks 
against local perturbations. If global perturbations are considered, situations can be dif- 
ferent. There is a possibility that the disk is stable against global perturbations. In the case 
of flickers in shell burning stars, the stars are unstable only against perturbations localized 
around the shell burning region (e.g., Schwarzschild & Harm 1965; Fujimoto, Sugimoto 
1979) and there is no instability for global perturbations. Non-negligible pressure change 
is associated with a global perturbation, and this acts to suppress any further temperature 
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increase in the shell burning region. Similar situations might be expected in geometrically 
thick disks if global perturbations are considered. 

Let us discuss the non-linear growth of perturbations in the case of a-viscosity. The 
growth rate of thermal instability is of the order of (aO)(A;r)"^/^ (if/r). The time interval 
during which a thermal perturbation propagates the radial distance r is rjU. Hence, the 
growth of perturbations during the time interval is ~ exp [{krfl^{Hlr)-^]. This suggests 
that the perturbations certainly grow during the propagation through the disk. A natural 
question is then, what happens at the non-linear stage. Do the thermal perturbations 
grow to a large amplitude or saturate with a small amplitude? To answer this question, 
the behavior of flickers in stars is suggestive. The flickers in shell-burning stars grow as 
long as the approximation of negligible pressure variation is valid. This approximation 
is violated flrst at a stage when the perturbation grows to a large amplitude and thus 
a pressure change in the shell-burning region occurs by a radial expansion of the whole 
star. After this the perturbation is damped out and the star returns to almost the same 
state as before. In other words, a sharp spike-like luminosity change occurs in the Kelvin 
time scale (i.e., thermal time scale). This sharp luminosity rise and fall is called a flicker. 
Flickers are repeated semi-regularly in a much longer time interval. We speculate that 
the local thermal instability discussed in this paper will behave like a flicker. That is, a 
local perturbation grows to a highly nonlinear stage where the pressure change becomes 
a non-negligible fraction. If so, we can expect spike-like time variations which will occur 
at various places of disks almost randomly in the thermal time scale. These perturbations 
propagate in the radial direction. This is interesting in relation to the observational facts 
that the rapid time variations of X-ray stars can be regarded as ensemble of shots (Negoro, 
Miyamoto, & Kitamoto 1994). 

The above arguments represent a picture of rather violent time variations. The a- 
model, however, may be irrelevant to represent the real situation, since the effects of 
a spatial variation of angular momentum distribution on the variation of stress tensor 
is represented only partially in this model. A diffusion-type viscosity will be better to 
represent the stress tensor due to turbulence. If we take a diffusion-type stress tensor, 
the growth rate of perturbations is of the order of aO as discussed in section 5.1, which is 
smaller than that in the case of the a-model. In this case, aO times r/?7 is about (if/r)~^, 
and perturbations may not grow much during their propagation through the disk. If so, 
the picture that the disks are covered randomly by small amplitude perturbations will be 
relevant. If the thermal diffusion process is considered, however, the growth becomes again 
rapid for optically thick advection-dominated disks. In optically thin advection-dominated 
disks local perturbations are damped. 
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Finally, we note limitations of our analyses and application of our results. The analyses 
in this paper have been made by using the vertically integrated equations. The use of such 
equations will be valid as the first approximation, if a) perturbations have no node in the 
vertical direction and b) the time scales associated with the perturbations are longer than 
the time scale of hydrostatic balance. These requirements result in a constraint on the 
viscosity parameter a, as discussed in section 2.1. To judge how much the approximation 
is really valid, however, a more careful consideration may be necessary, since we treat the 
case where the vertical scale of disks is much longer than the radial size of perturbations. 
To obtain a definite conclusion concerning this point, a comparison with two-dimensional 
analyses will be necessary. If a is not small, two-dimensional analyses will be necessary 
from the beginning in order to get reliable results. Stability against global perturbations 
is not examined in this paper. Furthermore, two temperature advection-dominated disks 
are also outside the scope of this paper. 

S. K. thanks Goteborg University and Chalmers University of Technology for their 
kind hospitality in 1994 May and June when this work was started. 



18 



References 



Abramowicz, M. A., Chen, X., Kato, S., Lasota, J. P., Ragev, O. 1995, ApJ, 438, L37 
Abramowicz, M. A., Chen, X., Krohk, J., Lasota, J. P., Narayan, R., Yi, I. 1995, in 
preparation 

Abramowicz, M. A., Czerny, B., Lasota, J. P., Szuszkiewicz, E. 1988, ApJ, 332, 646 
Abramowicz, M. A., Lasota, J. P., Xu, C. 1986, lAU Symp. 119, Quasars, ed. G. Swarup 

& V. K. Kapahi (Dordrecht: Reidel), 376 
Chen, X., Abramowicz, M. A., Lasota, J. P., Narayan, R., Yi, L 1995, ApJ, 443, L61 
Chen, X. M., Taam, R. 1993, ApJ, 412, 254 
Fujimoto, M. Y., Sugimoto, D. 1979, PASJ 31, 1 

Honma, F., Matsumoto, R., Kato, S., Abramowicz, M. A. 1991, PASJ, 43, 261 

Hoshi, R. 1977, Prog. Theor. Phys. 58, 1191 

Kato, S., Honma, F., Matsumoto, R. 1988, MNRAS, 231, 37 

Matsumoto, R., Kato, S., Fukue, J., Okazaki, A.T. 1984, PASJ, 36, 71 

Matsumoto, R., Kato, S., Fukue, J. 1985, in Theoretical Aspects on Structure, Activity, 

and Evolution of Galaxies III, ed. S. Aoki, M. lye, and Y. Yoshii (Tokyo Astronomical 

Obs: Tokyo), pi 02 
Narayan, R., Yi, L 1994, ApJ, 428, L13 
Narayan, R., Yi, L 1995a, ApJ, 444, 231 
Narayan, R., Yi, L 1995b, ApJ, in press. 
Narayan, R., Yi, L, Mahadevan, R. 1995, Nature, 374, 623 
Negoro, H., Miyamoto, S., Kitamoto, S. 1994, ApJ, 423, L127 
Pacyhski, B., Wiita, P. J. 1980, A&A, 88, 23 

Pedlosky, J. 1979, Geophysical Fluid Dynamics (Springer- Verlag, New York) 
Pringle, J. E. 1976, MNRAS, 177, 65 
Schwarzschild, M., Harm, R. 1965, ApJ, 142, 855 
Shakura, N. L, Sunyaev, R. A. 1973, A&A, 24, 337 
Wallinder, F. 1990, MNRAS, 237, 270 



19 



